Pressure-induced isostructural phase transition of metal-doped silicon clathrates 
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We propose an atomistic model for the pressure-induced isostructural phase transition of metal- 
doped silicon clathrates, Ba8Si46 and K8Si46, that has been observed at 14 GPa and 23 GPa, 
respectively. The model explains successfully the equation of state, transition pressure, change of 
Raman spectra and dependence on the doped cations as well as the effects of substituting Si(6c) 
atoms with noble metals. 
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Metal-doped silicon clathrates M8Si46 are the crys- 
tals that have the same structure (type I, space group 
Pm3n) as the methane hydrate (CH4)8(H20)46 [ij and 
the melanophlogite A2B6 (8102)45 where A and B are gas 
molecules 0. In the cubic unit cell, two small cages 
(Si2o) are composed of Si atoms at 16i and 24fc sites 
(see [1] for crystallographic notations); one cage is lo- 
cated at the corner and the other is at the body-center 
of the cell. The small cages are bridged by the direct 
covalent bonding between the Si(16i) atoms as well as 
by the covalent bonding via six additional Si atoms at 
the 6c sites, thereby forming six large cages (Si24) in the 
space between the small cages. Each small and large 
cage contains a metal atom M at its center, i.e., 2a 
sites and 6d sites, respectively. These metal-doped sil- 
icon clathrates are attracting much attention as nano- 
functional materials due to their remarkable properties 
in superconductivity 0, d, Hj Hi thermoelectricity d, Q, 
magnetism etc. The high-pressure study of these ma- 
terials are especially important because synthesis of va- 
rieties of these materials has become possible by using 
high pressure method 0, [ll| . 

Recently, a pressure-induced isostructural phase tran- 
sition has been reported for Ba8Si46 [H, [H, [l4] and for 
K8Si46 [H, [1^. On one hand. X-ray diffraction measure- 
ments show a sudden change of compressibility (soften- 
ing) in the lattice constant around 14 GPa for Ba8Si46 
[1^ and around 23 GPa for K8Si46 [III (see the symbols 
in Fig. [1]). Interestingly, the symmetry of the crystal is 
kept unchanged by this transition. On the other hand, 
Raman spectrum measurements show a drastic change in 
vibrational states around 15 GPa for B as 8146 (see Fig. 2 
of [1^) and around 20 GPa for K8Si46 [3, which prob- 
ably corresponds to the phase transition observed in the 
X-ray diffraction measurement. On the compression pro- 
cess of Ba8Si46 [iH the low energy peak about lOOcm^^ 
and some other peaks suddenly disappears at 15 GPa re- 



* Electronic address; ltiitaka@riken.jp] 

URL: http: //www. lltaka. orgT] 



suiting in six broadened Raman bands in high-pressure 
phase. On the decompression process the Raman spectra 
recover to those obtained before compression [l3| in spite 
of the drastic spectral change and the volume collapse. 
The disappearance and reappearance of the peak around 
100 cm~^ of Ba8Si46 at 15 GPa is an important signa- 
ture of this phase transition. In addition to this higher- 
pressure phase transition, another phase transition has 
been observed in Raman spectra at a lower pressure, i.e., 
7 GPa for Ba8Si46 J3\ and 15 GPa for K8Si46 flB^, whose 
signature seems absent in the X-ray diffraction [14]. The 
change of spectra at the lower-pressure phase transition 
is modest compared to the higher-pressure phase transi- 
tion, and is presently interpreted as the result of minor 
displacements of atoms from their ideal position [3, [13] ■ 
Mechanism for the higher-pressure phase transition have 
been proposed as a change of hybridization between Ba 
atoms and the silicon cages by San Miguel et al. and 
as an electronic topological transition (ETT) by J.S. Tse 
et al. [l3|, however, the observed collapse of the lattice 
parameter has not yet been explained quantitatively by 
these mechanisms. 

In this article, we propose a new model that explains 
many aspects of the observed data such as lattice con- 
stants, Raman spectra, and the effect of noble metal sub- 
stitution. Our model assumes that vacancies start form- 
ing in the silicon framework, say at Si(6c) sites, above the 
transition pressure ( e.g. 15 GPa for Ba). Si(6c) sites are 
partially occupied by silicon atoms with random config- 
urations, and the crystal transforms to Ba8Si46-n, where 
< n < 6 and typically n ~ 3 at the highest pressures, 
i.e., 

Ba8Si46 Ba8Si46-n + "-Si. (1) 

The assumption of random configuration implies that, on 
the average, only the occupation number of Si(6c) sites 
changes to (6 — n)/6 and the space group Pm3n is pre- 
served. On one hand, however, the crystal may become 
more compressible due to the created vacancies, result- 
ing in the collapse of the lattice parameter. On the other 
hand, the partial occupancy of Si(6c) sites implies weaker 
binding between small cages, which may result in the 
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FIG. 1; (a) Equation of state for the nominal BagSi46. Sym- 
bols indicate the experimental data [Ij. Solid lines indicate 
the theoretical results for Ba8Si46, and the dashed lines for 
Ba8Si43. (b) Equation of state for the nominal K8Si46. Sym- 
bols indicate the experimental data 0]. Solid lines indicate 
the theoretical results for KgSi46, and the dashed lines for 

K8Si45. 



softening of some Raman modes. These features will be 
confirmed in the following by the first principles calcula- 
tions. Actually, this type of partial occupation has been 
observed for the nominal Ba8Ge46, an analogue material 
of Ba8Si46, at ambient conditions The favorable 

stoichiometry is reported as Ba8Ge43 |18j . 

To examine the validity of our model, we first made 
geometric optimization of the unit cell for Ba8Si46-n, 
n = 0,1,3,6 under pressure. The first-principles calcu- 
lations were performed with CASTEP 4.2 codes based 
on the plane wave basis set[l^, the Vanderbilt-type ul- 
trasoft pseudopotentials [l^] for electron-ion interaction, 
and GGA-PBE ^21j for exchange-correlation interaction. 
The cut off energy was set to 160 (eV) and the Brillouin 
zones were sampled with 2x2x2 k-points. For n = 1 
(n = 3), the partially occupied Si (6c) sites with random 



FIG. 2: (a) Enthalpy of BagSi4o, Ba8Si43 and Ba8Si45 rela- 
tive to Ba8Si46 estimated by first principles calculation, (b) 
Enthalpy of K8Si43, K8Si44 and K8Si45 relative to K8Si46- 



configurations were further modeled by the ordered unit 
cell with one (three) out of six 6c sites are occupied by 
Si atoms with maximum symmetry. As the result, the 
unit cell for n = 1 becomes tetragonal with PAm2 sym- 
metry, and for n — ?> trigonal with i?32 symmetry. After 
geometric optimization, it was found that these unit cell 
is close to the original PmSn symmetry and may ap- 
proximate random occupation well. Symbols in Fig. [1^ 
show the measured equation of state (EOS) of the nom- 
inal Ba8Si46 [ij]- The measured EOS starts to deviate 
from the calculated EOS of Ba8Si46 around 15 GPa and 
approaches to the calculated EOS of Ba8Si43 above 20 
GPa. This result strongly supports our scenario men- 
tioned above. 

Next, we study the enthalpy for vacancy formation to 
understand why vacancies start forming around 15 GPa 
and why n approaches to ~ 3. For this purpose, we define 
the enthalpy of Ba8Si4o, Ba8Si43 and Ba8Si45 relative to 
Ba8Si46 as follows. 

A7?(Ba8Si46-„) = F(Ba8Si46-„)+ nH {Si) -HiBasShe) 

(2) 



3 



where i? (Ba8Si46-n) is the enthalpy of Ba8Si46-ri, and 
iJ(Si) is the enthalpy per atom of the diamond silicon 
(below 10 GPa) and silicon V phase (above 10 GPa) 
[23 , [l^l. For simplicity, we have ignored the interme- 
diate minor phases between the two phases. The results 
are shown in Fig. [2|i. AiJ(Ba8Si45) becomes negative 
around 15 GPa, indicating the onset of vacancy forma- 
tion at Si(6c) sites. Then A_ff(Ba8Si43) becomes nega- 
tive around 22 GPa, while Ai?(Ba8Si4o) remains posi- 
tive even above 30 GPa. Therefore vacancies start form- 
ing around 15 GPa, and then Si(6c) sites approaches to 
the half occupation but will not be fully unoccupied. 
When the pressure changes, the number of vacancies 
at Si(6c) sites is adjusted to minimumize the enthalpy 
Ai7(Ba8Si46-n) • This scenario explains the phase tran- 
sition of Ba8Si46 measured at 15 GPa |1J] very well. This 
phase transition may be induced by the phase transition 
of the material composing the clathrate, i.e.. Si-diamond 
— > Si-V; Si atoms are squeezed from the clathrate be- 
cause the more stable Si phase appeared. Similar phe- 
nomenon is known for the pressure-induced phase transi- 
tion of methane hydrate '24, 2^ [H, [13, [111 where water 
molecules are squeezed from the hydrate crystal in ac- 
cordance with the pressure-induced phase transitions of 
water ice. 

The electronic configuration of atomic Ba is 
[Kr](4d)i°(5s)2(5p)6(5d)"(6s)2. However, Ba 6s states 
are not occupied in Ba8Si46: part of these two electrons 
occupy the hybridized states of Ba 5d and Si 3p states 
that form the narrow DOS peak at the Fermi energy, 
and the rest of the electrons are donated to the Si 
framework p9| . Therefore Ba is partially ionized. When 
one Si atom is removed from the crystal by forming 
a vacancy, four dangling bonds are formed on the Si 
framework, and these dangling bonds should be satisfied 
by consuming four electrons transferred from the Ba 
atoms [l^. Since there are eight Ba atoms in the unit 
cell and each Ba atom can donate two '6s' electrons, 
the maximum number of vacancy formation per unit 
cell becomes (8 x 2)/4 — 4. Creating more vacancies 
will leave unsaturated dangling bonds, resulting in 
much higher energy of formation. The extremely high 
enthalpy of Ba8Si4o compared to Ba8Si45 and Ba8Si43 
in Fig. [2^ can be explained by this mechanism. Further, 
the increase of the width and the area of Ba Lm X-ray 
absorption near-edge white line spectra as a function 
of pressure above 11 GPa 44] may be interpreted as 
a result of vacancies randomly generated on the Si 
framework. 

Now, let us discuss vacancy formation of K8Si46. Sym- 
bols in Fig. [TJs show the measured EOS of the nominal 
K8Si46 .9]. The measured EOS starts to deviate from the 
calculated EOS of K8Si46 around 20 GPa and approaches 
to the calculated EOS of K8Si45 around 30 GPa. Re- 
member that the electronic configuration of atomic K is 
[Ar](4s)"'^. Then the maximum number of vacancy for- 
mation per unit cell becomes (8 x l)/4 — 2. Fig. [^h 
shows the enthalpy AiJ(K8Si46-n) as a function of pres- 



sure. It predicts that vacancy formation starts at 20 GPa 
and the crystal remains approximately K8Si45 up to 30 
GPa. This result is consistent with the measured EOS. 
The degree of volume collapse seems controlled by va- 
lence number of the cation. 

In order to explain the observed Raman spectra, we 
calculated the normal modes of Ba8Si46 by using density 
functional linear response theory [s^, [3l|, [3^ . For each 
mode, we examined the symmetric characters and the 
contributions from each atoms. The mode around 100 
cm~^, which is an important signature of the phase tran- 
sition at 15 GPa, consists of small cages rotating around 
their center, which are then tied up to a collective mode 
being bridged by Si(6c) and Si(16i) atoms. According to 
our assumption. Si atoms are partially and randomly oc- 
cupied above 15 GPa, and therefore the restoration force 
of this collective mode will become weaker and randomly 
distributed, resulting in a red shifted and broadened 
peak. The observed disappearance of the peak around 
100 cm^^ and appearance of the broad peak around 60 
cm~^ above 15 GPa may be interpreted as the result of 
the partial occupation of the bridging atoms. 

Ternary Si clathrate Ba8Au6Si4o is synthesized under 
high pressure from the stoichiometric mixture of Ba, Au 
and Si, where Si(6c) atoms of Ba8Si46 are replaced by 
Au atoms [l^. This indicates that the Si (6c) atoms are 
weakly bound to the framework, and replacing Si(6c) 
by Au atoms makes the system more stable. Our first 
principles calculations of enthalpy show that Ba8Au6Si4o 
is 5 '-^ 10 eV more stable than Ba8Si46 in this pres- 
sure range. According to our model, this stabilization of 
6c atoms should suppress the creation of vacancies and 
hence the isostructural phase transition. Actually, recent 
Raman spectra measurement revealed that Ba8Au6Si4o 
is stable up to 27 GPa with no phase transition ob- 
served 33]. Similar stabilization of Ba8Ag6Si4o has been 
observed by X-ray diffraction (33 |. 

So far we have discussed the energetics of the isostruc- 
tural transformation. Let us now examine the kinetics of 
this transformation. According to our model the excess 
Si atoms are transported by vacancy mechanism. Va- 
cancies are created at the surface or grain boundary of 
Ba8Si46 particles by releasing the excess Si atoms from 
the particle and then the created vacancies propagate 
into the particle. In Fig. [31 the enthalpy of a vacancy 
at 15 GPa is shown as a function of position along the 
covalent framework, where V(6c) is set to the origin of 
enthalpy. In the graph, for example, the vacancy moving 
from the 6c site to the 24fc site means the Si atom at 24fc 
site moves to the empty 6c site. Since equation Q defines 
A7?(Ba8Si45) = V{6c), vacancies start to form at site r 
when Ai7(Ba8Si45) + V{r) < 0. In contrast to Ba8Ge43, 
it turned out that vacancies favor energetically the 16i 
and 24k sites than 6c sites. When vacancies start to form 
at these crystallographic sites, the migration mechanism 
of vacancies is the thermal hopping between the sites sep- 
arated by the potential barriers of c.a. 0.2 eV, which is 
smaller than the 0.4 eV barrier of the diamond silicon at 
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FIG. 3: Enthalpy of a vacancy relative to the formation en- 
thalpy at 6c site ( P—15 GPa ) is presented along the covalent 
framework. 

ambient pressure [s^ but still large enough to make the 
diffusion slow at ambient temperature. However, once 
the enthalpy difference becomes large enough, 

Ai7(Ba8Si45) < -O.l(eV), (3) 

to override the barriers, vacancies can propagate ballis- 
tic into the particle. This simple mechanism may explain 
the strong mobility and diffusion necessary for our model, 
while more realistic models may take into account the ef- 
fects of energy dissipation and many vacancies. Since the 



condition ([3]) is almost equivalent to the vacancy forma- 
tion condition for 6c sites, Ai?(BagSi45) < 0, the discus- 
sions made with the assumption that 6c site is the most 
favorable vacancy position, are still valid in this case. 

In summary, by using the first principles calculations, 
we have shown that the model of partially occupied Si 
sites can explain many aspects of the phase transition of 
Ba8Si46 at 15 GPa, such as the equation of state, tran- 
sition pressure and change of Raman spectra as well as 
effects of substituting Si(6c) atoms with noble metals. 
Recently, Carrillo-Cabrera et al. [s^ and Okamoto et 
al. [13] have reported ordered models of Ba8Ge43 with 
2x2x2 superstructure. Our discussions of vacancy for- 
mation should be valid also for such superstructures. For 
further confirming, improving, or discarding our model, 
more experimental and theoretical studies are necessary, 
especially on the reversibility of the phase transition and 
the fate of the released Si atoms. Precise X-ray diffrac- 
tion measurement may provide a direct evidence of the 
partially occupied Si sites or that of the superstructures. 
Raman spectrum measurement of the analogue material 
BagGe43 at ambient and high pressure conditions will 
provide a clue to interpret the Raman spectra of the nom- 
inal Ba8Si46 above 15 GPa as those of Ba8Si43 phase. 
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